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Abstract 

We have investigated KN-ttY system with a coupled-channel complex scaling 
method (ccCSM). By an advanced use of the ccCSM, scattering states as well 
as resonant states in the system have been treated in the same framework. 
Respecting the KN scattering length obtained by the Martin's analysis, we 
have constructed a meson-baryon potential matrix based on a chiral SU(3) 
theory, with a local Gaussian form in r-space. To test the property of the 
constructed potential, we have studied the scattering amplitude for isospin 
1 = and 1 = 1 channels and a resonant state in the 1 = channel, in the 
ccCSM with both non- and semi-relativistic kinematics. The resonant pole 
corresponding to A(1405) is found around (1419, -20) MeV ((1419, -13) 
MeV) on complex-energy plane with the non-relativistic (semi-relativistic) 
kinematics. Mean distance between meson and baryon in the resonant state 
is about 1.3 — i0.3 (1.2 — z0.5) fm, which means that the state is considered 
as a Gamow state. We have observed a signature of another pole around 
(~ 1360, —40 ~ —90) MeV, although it is unstable against change of scaling 
angle 9. This may correspond to the lower pole of the double-pole of A(1405) 
discussed in literature today. 

Keywords: KN-ttY system, A(1405), scattering amplitude, complex 
scaling method, chiral SU(3) theory 
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1. Introduction 



^-nuclear system has been a hot topic in nuclear and hadron physics for 
a long time. Due to strongly attractive KN interaction in isospin 1 = 
channel, finite nuclear systems with anti-kaons are expected to have exotic 
properties such as deeply bound and quasi-stable states with high density 
|2J. Such kaonic nuclei have been investigated with various many-body 
treatments 0, 0, H|. In particular, to clarify the property of kaonic nuclei, 
great efforts have been devoted to investigate f^pjfl, a prototype of kaonic 
nuclei, in both of theoretical and experimental studies. From the theoretical 
studies, it have been claimed that K~pp will not be so deeply bound and 
its decay width will be large, (total binding energy < 100 MeV and decay 



width > 50 MeV) 
between calculations 



although there are quantitative discrepancies 
On the other hand, experimental results indicate 



deeper binding of K pp than theoretical predictions if the observed state is 



PJL 



111 ], although there are some objections to the exper- 
Thus, the consensus for K~pp has not been achieved 



the bound K~ 
imental result 
yet. 

For study of kaonic nuclei, KN (involving tcY) interaction is a basic 
input, and a resonant A(1405) is an essential building block because it can 
be reasonably interpreted as a quasi-bound 1 = state of KN with s-wave, 
rather than a three-quark state [l3|. One approach to KN system is a study 
based on the chiral SU(3) coupled-channel dynamics [lij]. This approach 
(called the chiral unitary model by Oset and Ramos 15j]) has been succeeded 
in studies of s-wave meson-baryon systems including S = —1 sector 16( . The 
double pole nature of A(1405) pointed out within this model, is interesting 
171 . 1 1 81 ] . According to further studies along this model, experimental data 
seem to support the double-pole nature [l9|. Recently, accurate data at the 



KN threshold are given by a precise measurement of Is level shift of kaonic 



hydrogen atom [20 
region 



in addition to experimental data on KN subthreshold 



21(. Due to such precise data, the physical quantities near the KN 



threshold are strictly constrained and those uncertainties below the KN 



threshold are also expected to be decreased [22 . 

In such a current situation, we start a study of kaonic nuclei with a 
coupled-channel complex scaling method (ccCSM), keeping in mind the fol- 



1 Actually, this system is a KNN-ttYN coupled system with quantum numbers J 71 " = 
CP, (T,T Z ) — (1/2, 1/2). It is expressed symbolically as K~pp. 
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lowing three points: 1. Simple and adequate treatment of resonant states in 
many-body system, 2. Explicit inclusion of all channels in a coupled- channel 
problem, and 3. Accessible to structure of kaonic nuclei. Looking back past 
studies on K~pp, in variational studies [H, 0] it is treated as a bound state 
in KNN channel as a consequence of the elimination of ttY channels. In 
Faddeev-AGS studies (HI, @] certainly a coupled-channel calculation is fully 
performed and resonant poles are searched. However, there a separable form 
is assumed to the KN-ttY potential which is a key ingredient in the K nu- 
clear study, and a wave function is not obtained explicitly in this approach 
though it is important to investigate the nature of kaonic nuclei. Thus, each 
approach involves advantages and disadvantages. Since the ccCSM is ex- 
pected to overcome above disadvantages, we employ this method to study 
the K nuclear system. 

The complex scaling method (CSM) has been applied to various nuclear 
physics, and it has been greatly succeeded in particular in the study of res- 
onant states of unstable nuclei jiij. The CSM is a practical tool for the 
study of nuclear many-body systems. Indeed, a resonant nature in unsta- 



ble He is revealed with the CSM in Ref. [24| where a five-body system of 
4 He + n + n + n + n is solved. The CSM is suitable for resonant states, 
because we can handle them in the same way as bound states. Though a 
resonant wave function is originally divergent at infinite distance, in the CSM 
it is transformed to a square-integrable function by a complex rotation for 
the coordinate and then it can be represented with e.g. a Gaussian base 
which is familiar for ones study bound states. In addition, by an advanced 
used of complex-rotated wave functions, the scattering amplitude can also 



be calculated with the Gaussian base [25(. Thus, all of bound, resonant 
and scattering states can be handled in a single framework of the CSM with 
Gaussian base. 

Since it is our first attempt applying the ccCSM to ^-nuclear systems, 
in this article we investigate the s-wave two-body system of KN-7rY coupled 
channels. We examine semi-relativistic kinematics as well as non-relativistic 
one to be careful of a pion which is a light-mass particle. First, we will 
check how the ccCSM works in the present system. Then, we will construct 
a meson-baryon potential for the KN-ttY coupled system, based on a chiral 
SU(3) theory. We adopt a Gaussian form in the coordinate space to the 
potential shape, because of convenience for our further study of ^-nuclear 
system with Gaussian base. Our potential is constrained by the KN scatter- 
ing lengths for both isospin states obtained by the Martin's analysis of old 
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KN scattering data [26|. Using the constructed potential, we will investigate 
the behavior of the scattering amplitude. In the 1 = sector, the pole on the 
complex-energy plane also will be investigated in detail because there should 
be a resonance corresponding to the A (1405). 

The contents of the article are as follows: In the section 2, we will explain 
a meson-baryon potential used in our study and the formalism of ccCSM 
for the study of scattering state as well as resonant state in detail. In the 
section 3, the obtained results will be given. The scattering amplitudes will 
be shown for both isospin states. For the 1 = channel, the structure of a 
resonant state will be investigated. In the last section, we will mention our 
summary and future plans including some discussions. 

2. Formalism 

2.1. Kinematics and interaction 

We are considering a KN-nY coupled system in s-wave, where Y in- 
dicates a hyperon which is S (A and S) for / = (J = 1) case. We in- 
vestigate such a two-body system in semi-relativistic kinematics as well as 
non-relativistic one, because a pion joins in our calculation and its mass is 
very small. The Hamiltonian for non-relativistic kinematics is 

H = ^[M a + m a + ^ |a)(a| + V$§, (1) 

and that for semi-relativistic kinematics is 

{^/M^+¥ a + + \<*){a\ + Vg&. (2) 

a 

Here M a , m a and baryon, meson and a reduced mass in the channel 

a, respectively. p a is the relative momentum between a meson and a baryon 
in the channel a. 

The last term V^§ {Ymb) represents a meson-baryon potential for non- 



relativistic (semi-relativistic) kinematics. In Ref. [14] . Kaiser, Siegel and 
Weise proposed a pseudo-potential for meson-baryon system with S = — 1 
derived from an effective chiral SU(3) Lagrangian, and used it in a Lippmann- 
Schwinger equation to obtain A(1405). In this paper we follow this work and 
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adopt the lowest order of the pseudo-potential: 



Vmb = E-^ u ° + u e)\^&{r)\aW\, (3) 

Stp(r) = (v^^)" 3 exp[-(r/^) 2 ], (4) 

where the spatial dependence is modified from the 5 function in the origi- 
nal to a Gaussian function with a range parameter d a g. Hereafter, we de- 
note this as "KSW potential". Note that this potential is energy depen- 
dent because meson energy u a , baryon energy E a and the reduced energy 
Cj a = u a E a /(u a + E a ) are given by function of CM energy y/s as follows: 

2y/S 2y/S 

Note also that structure of this potential is determined by the coefficients 
{Cap} which are computed by Clebsch-Gordan coefficients of SU(3) and given 
as 



KN ttE 




7rA 



c d=o) = ( * "V )» ^ (/=1) = I 2 J. (6) 

Strength of this potential depends on the pion decay constant /„-. It is noted 
that the range parameter d^g can be a different value in each channel set 
(at,/?). For the range parameter of the transition potential, we assume that 
d a p = (d aa + dl g )/2 to reduce the number of parameters. 

Since the flux factor ,/ MaA fe in Eq. (EI) is based on the relativistic kine- 

matics, the KSW potential had better be used in the semi-relativistic frame- 
work in our calculation. When we adopt the non-relativistic kinematics, 
we make non-relativistic reduction to the flux factor in two prescriptions, 
denoted as "NRvl" and "NRv2". In the first prescription, comparing the 
expression of the differential cross section in both kinematics, we replace the 
reduced energy in Eq. (J3j) with the reduced mass; 



NRvl: V^r = Y,-%^ + ^\^3Ur)\a)^\. (7) 
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In the second prescription, by considering the small momentum limit, we 
replace the meson and baryon energies with the meson and baryon masses, 
respectively; 



NRv2: C 2 = E-%( w « + ^ltVi«H(^ W 

a,/3 J7T V a ' & 

In both non-relativistic approximations, we keep remaining the meson-energy 
part (w a + ojp) as original, because this energy dependence is attributed to 
the chiral dynamics which we respect in our study. Of course, energies such 
as a/s and u a themselves, are calculated by the non-relativistic formula. 

2.2. Coupled- channel complex scaling method 

As mentioned in the introduction, the s-wave KN-tcTj system has a res- 
onant state in / = channel which corresponds to the A(1405). Such a 
resonant state of a meson-baryon system can also be investigated with the 
complex scaling method (CSM), in the same way as resonant states of unsta- 
ble nuclei. Here, we give a brief explanation of the coupled-channel complex 
scaling method (ccCSM). Details of the CSM and its successful application 



to unstable nuclear physics are summarized in Ref. [23 

In the CSM, a relative coordinate r and a conjugate wave number k in 
Hamiltonian H and a wave function |$) are complex-scaled as 

U{6) : r -»■ re ie , k -»■ fee - * (9) 

Then, transformed Hamiltonian and wave function are defined as He = 
U(6)HU^ 1 (9) and |$6») = U(8)\$), respectively. We expand the wave func- 
tion <3>(r) in partial waves as $(r) = J2im 4 ) i( r )/ r ^zm(^) as usual. In our 
definition, all radial wave functions are transformed by the complex scaling 
operator as 



#(r) = U{e)Ur) = e^ 2 Mre t0 ), (10) 

taking into account the Jacobian in integration to calculate expectation value 
of operators. An expectation value of operator O is calculated with bi- 
orthogonal set {&e}: 

(6)e = ($e\de\Qe), (11) 
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where $g(r) = <&g{r) in terms of radial part of wave function for bound and 
resonant states. The complex scaled wave function is generally normalized 
as (&e\&g) = 1. With the definition of the bi-orthogonal state, the radial 
part of the complex scaled wave function is normalized as 

POO 

£ / dr{#(r)} 2 = l (12) 
i Jo 

for bound and resonant states [23J. 

It is known that energies of bound and resonant states are independent 
of the scaling angle 9 while those of continuum states vary with 9 as ^e 2ld 
in case of non-relativistic kinematics. (ABC theorem [23[) In addition, it 
is easily understood that a resonant wave function is transformed from a 
divergent function to a dumping function by the complex scaling U(9) with 
adequate values of 9. In other words, the boundary condition for resonant 
states is modified to the same one for bound states. Therefore, we can obtain 
resonant states as follows: For various #'s, complex eigenvalues are calculated 
by diagonalizing the complex-scaled Hamiltonian Hq with Gaussian base as 
done in usual studies of bound states. Among obtained eigenstates, the 
states with the eigenvalues independent of 9 are recognized as resonant states. 
Since the continuum states appear along a line tan~ 1 (Ei/En) = —29 on the 
complex-energy plane, resonant states can be separated from the continuum 
states when we set the scaling angle 9 appropriately. Similarly, resonant 
states can be found also in case of semi-relativistic kinematics, although 
complex eigenvalues of continuum states have different ^-dependence from 
non-relativistic case. 

In the present study, the CSM is applied to a coupled-channel problem, 
because the 1 = wave function contains the KN and 7rS components. The 
spatial parts of these wave functions with s-wave, 0/nv(r) an d 0vrs(^), are 
expanded with Gaussian base G%(r): 

\^rn- = 1) = l<f>RN(r)Yoo(mKN) + ±<fi*(r)Y M {n)\TrV) (13) 



kl=0i 



]Tq^(r) = J2 C f N l= ,{b 3 )reM-r 2 /2b% (14) 



where coefficients {C°} are complex parameters to be determined. As ex- 
plained above, by diagonalizing the Hg with the basis {G"(r)\a)} and in- 
vestigating the 9 dependence of eigenvalues, we can find resonant states. As 
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explained in Eq. (TT2"j) . the radial part of the 1 = wave function (Eq. (TT5T) ) 
is normalized as 



do 

ill 







(15) 



when it is complex-scaled. 

Note that the self-consistency for the complex energy is needed to be 
considered, since the meson-baryon potentials (Eqs. (J3J), (j7|) and (jSJ)) have 
an energy dependence attributed to chiral dynamics. The eigen energy E ca i c 
of Hamiltonian calculated with the ccCSM should coincide with the energy 
E int inputted into the meson-baryon potential as Vmb(Vs = E int ). Here 
E ca i c and Ei nt are complex values because they are the energies of resonant 
states, Er — iY/2. (Er and T are energy and decay width, respectively) 
At the n-th iteration, we set a value of E^j as an input for the Vmb and 
then obtain an eigen energy E^j c of resonant state with ccCSM. We use the 
E^j c as an input E^ 1 ^ for the next tern. Such iterations are repeated until 
the self-consistent condition is satisfied; E^j c = E^j. In case of the present 



system, the self-consistency is achieved in five-times iterations 27|. (n < 5) 



2.3. Scattering amplitude calculated with ccCSM — "CS-WF" method— 

We are interested in the scattering state as well as the resonant state 
which can be investigated with the usual CSM as explained in the previous 
section. In this article, we solve a scattering problem also with a square- 
integrable base such as Gaussian base with the help of the CSM, following a 
method called "CS-WF" which was developed by Kruppa, Suzuki and Kato 



251 ] . Here, the detailed formalism of CS-WF is shown for a multi-channel 



case such as the KN-ttY system. 

We assume a system involving n channels in which the channel Co is the 
incident channel with the energy E. In a non-relativistic case, the radial 
Schrodinger equation for such a multi-channel system is 

n 

(E - H l c )^\r) = £ V^ir), (16) 

c'=l 

, h 2 d 2 h 2 i(i + i) Tr , , „, 

where \i c and Mr )C are the reduced mass and the total mass in the channel 
c, respectively. c (r) is a direct potential for the channel c and V cc i is 
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a transition potential between channels c and d . We assign closed (open) 
channels for the incident energy E to channel numbers c=l~n^(c = 
tib + 1 ~ n). The wave function of the channel c can be written as 

,W W J 4f' B (r) (c=l,---,n B ) 
l ' c{> 1 Uk Co r) 5 C , C0 + 4f s \r) (c = n B + 1, • • • , n) ' 1 1 

The incident wave number k rn satisfies 7r—k^ + Mt„. = E. 

Inserting this equation into Eq. ffToT) and applying the complex scaling, 
then we obtain the coupled equations as follows: 



(E - Hjf) 



1 ( r ) 



(c ),sc,9 



r 



" B 



c'=l c'=n B + l 



where ^c°' ) ' S ' 6 '( r ) (W ^°- ) ' sc ' 6 '(r)) means the i^i^' B (r) (V"jc ,sc ( r )) complex-scaled 
as following Eq. ( ITU]) . The complex-scaled Hamiltonian for the channel c is 



^ = -e-™^-^- 2 + e -^pl+A + ^ (re «) + M (20) 
2ji c ar l 2fi c r 2 



The wave functions for the closed channels, ^ S (r)}, are square- 
integrable since the threshold energy of these channels is above the incident 
energy E. The complex-scaled ones are also square-integrable. As for the 
open channels whose threshold energies are below the E, the scattered part 
of wave functions, ^ sc (r)}, are not square-integrable, since they behave 
in the asymptotic region as 

W,c 0) ' S >) -> k c fi,cc (kc)H +) (k c r) oc exp{*(A; c r - lix/2)} at r -> oo, (21) 

where hf(kr) is the Riccati-Hankel function. But, they are transformed to 
be square-integrable functions due to the complex scaling. The asymptotic 
behavior of the complex-scaled scattered part of wave function is 

ipi^' sc ' e (r) oc i~ l exp {ik c r cos 6 — k c r sin 6} at r — > oo. (22) 



It is easy to understand that the ip^'' sc ' e becomes a square-integrable func- 
tion for < 9 < 7i. Thus, since both the complex-scaled wave functions 
{Wc°' > ' B ' 6, ( r )} an d {Wc°' ) ' SC ' 6 '( r )} are square-integrable, they can be expanded 
with Gaussian base {Gj(r)} as 



N 



(23) 



In the present study, we use a common set of normalized Gaussians for all 
channels: 



Gi(r) = N l {b l )r l+l exp 



2b? 



Nfc) = 6r (aw,)/a 



2 Z+2 



(2/ + 1)!! 



• (24) 



Inserting the Eq. (123]) into the coupled equations f lT9|) . linear equations 
for the unknown variables {t^f(9)} are obtained: 



J L 



CCQ,l 1 



(25) 



where each matrix element indicates Oy = = \Gj), 

= mV^re^lGj) and " 



6° • = e" 



^ 2 / drG,(r)y cco (re^)j z (^). 



(26) 



We explain how to calculate scattering amplitudes in the remaining part 
of this section. With the scattering wave functions {^\ c °\ r )}j the scattering 
amplitude between the initial channel cq and the final channel c is expressed 
as 



fl,cc {kc) 



fJ>c 



h 2 k c k C0 c , =1 jo 



n poo 

£ / drUk c r)V cd {r)^{r). (27) 



It is obtained with the help of Green function, as a detailed explanation 
Appendix A. By inserting Eq. ffTSj) into the $ j'(r) of the 



is given m 



above equation, the full scattering amplitude is decomposed to the Born 
term /j^ n (& c ) attributed to the incoming wave ji(k Co r) and the other part 
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ficcoi^c) attributed to the scattered wave ip\ c ° (r). The Born term can be 
obtained by the numerical integration. In the calculation of fi S c Co {kc), the 
integration path can be modified from r-axis to re^-line due to the Cauchy's 
theorem. Therefore, the /f£c (& c ) is equal to f[ C c 'c (,k c ) that can be calculated 

with the complex-scaled wave functions (r)}, as 

. 71 poo 

&c) = -J£- e' d/2 £ / * k**e*) V c Are*) v£°*(r). (28) 
a tv c tv co Jo 

Thus, we can obtain the full scattering amplitudes as 

Accote) = Cc?(k c ) + fr c c (k c ) (29) 
fl B Z n (kc) = -J^t- r drUk c r)V cco {r)Ukc r) (30) 

n K C K C0 Jo 

n N 

fi s Uk c ) = /Sw^-^rEE^wt- ( 31 ) 

C C c ' = l j=l 

The filf(k c ) can be obtained with the matrix elements {b e c , ci } by using the 
Eqs. d#D and $2B). 

We summarize the essential points of the CS-WF at the end of this sec- 
tion. The first point is that the incoming part ji(k co r) is separated from the 
scattered wave functions {^j c ( r )} as snown m Eq. ( !T8|) and that only the 
scattered parts {ip[ C c\ r )} are complex-scaled. The complex scaling is used 
to make a non-square integrable function transformed to a square integrable 
one. If the full scattering wave function $^ \r) is complex-scaled, it does 
not become a square-integrable function because of Riccati-Bessel function 
ji(k CQ r) that contains both components of exp{±i(/cr — ^)} m the asymptotic 
region. The second point is that we calculate the scattering amplitude with 
the complex-scaled function ip^ ( r )> instead of the ipi° (r) which is needed 
in usual calculation of the scattering amplitude. By the Cauchy's theorem, 
the amplitude can be obtained with the if/j^ (r) which is expressed with a 
square-integrable base. We note that the scattering amplitudes calculated in 
this way are independent of the scaling angle 9. 

In case of the semi-relativistic kinematics, the scattering amplitudes are 
calculated in the same way as the non-relativistic kinematics as explained 
above, by replacing the matrix elements of kinetic energy term with those 
for the semi-relativistic kinematics. These matrix elements are given in 
Appendix B 
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KbarN 



E Kbar N [ MeV ] 



Figure 1: (Left) Unitarity violation of the S'-matrix. (Right) Comparison of the phase-shift 
sum between CS-WF (red diamond) and CLD (black line). "ExbarN" means the energy 
of the system measured from the KN threshold. Here, the 1 = channel is considered. 



2.4- Test of ccCSM for scattering amplitude 

We test the CS-WF method in I = KN-tcY, system, since it is the 
first time to apply this method to a meson-baryon system. Here, we employ 
a phenomenological KN potential [l[ which is a coupled KN-nT, potential 
and a local potential with a single Gaussian form in r-space, but an energy- 
independent potential. 

First, we check the unitarity violation of S'-matrix; 1 1 det 5" | — 1|. We 
calculate the S'-matrix for the partial wave I with the scattering amplitude 
obtained by the CS-WF method. The S-matrix element for channels a and 
(3 is related to the scattering amplitude as 



SL = ^ + ^J^f^ a , (32) 



where k a is the wave number in the channel a. As shown in the left panel 
of Fig. [H the magnitude of the unitarity violation is confirmed to be sig- 
nificantly small of the order of 10~ 5 in KN energy region of —120 to 30 
MeV. 

We check also the phase-shift sum. We compare the sum of phase shift, 
^kat + ^ttE, with that obtained by the continuum level density (CLD) method 



281 ] . Noted that the CLD method can give the phase-shift sum by using 
eigenvalues of the complex-scaled Hamiltonian, but that it cannot give the 
phase shift of each channel. Phase shift is obtained from the S-matrix (Eq. 
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-150 -100 -50 50 

E K barN [ MeV ] 

Figure 2: The KN — > KN scattering amplitude in the I — channel calculated with a 
phcnomcnological potential (Jj . 



(PD) as 

q= fS n S 12 \ = f cos2ee i51 % sin e e l ^ +5 ^ \ , , 

\S 21 S 22 J \ ismee^ 1 ^ cos2ee i52 )' 1 ' 

where 5i is phase shift of channel i and e is a mixing parameter. As seen 
in the right panel of Fig. (TJ the phase-shift sums calculated with both the 
methods agree with each other quite well. 

We show the KN scattering amplitude calculated with our method in Fig. 
|2j The KN scattering amplitude is given in Ref. [l[ where it is obtained by 
solving the Schrodinger equation as usual. Compared them with each other, 
we find that our calculation reproduces well the original result. (See Fig. 
1 in Ref. The scattering length also coincides well in both calculations: 
— 1.77 + i0.47 fm in our calculation and — 1.76 + i0.46 fm in the original work. 



3. Result 

We show the results of our calculation with the ccCSM. As mentioned 
in the section 12.11 the KSW potential and two kinds of non-relativistically 
approximated one are used, and two kinematics, non-relativistic and semi- 
relativistic, are examined. Hereafter we abbreviate non-relativistic and semi- 
relativistic as "NR" and "SR", respectively. We consider natural combina- 
tions of potential types and kinematics as follows: 
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- Case (i) One of non-relativistic approximated KSW potentials 

(Eq. (j7])), NRvl, with the non-relativistic kinematics. 

- Case (ii) The other non-relativistic approximated KSW potentials 

(Eq. OH])), NRv2, with the non-relativistic kinematics. 

- Case (iii) The KSW potential (Eq. 03])) with the semi-relativistic 

kinematics. 

In addition, we investigate an instructive case that the KSW potential is 
used with the non-relativistic kinematics. 

We consider f n in the KSW potential as a parameter in our model. In 
the present study, it is varied around the physical values of f n ~ 93 MeV and 
f K ~ 110 MeV. We examine four cases of f w = 90, 100, 110 and 120 MeV. 

3.1. Scattering amplitude of I = KN-nY, system 

First, we determine the range parameters of Gaussian functions, {d 1 ^ }, 
in the meson-baryon potentials, Eqs. ([3]), ([7]) and (]8]). In the present study, 
with an assumption of d^N,^ — d^.RN — (^i?7v,i?7v + ^7rS,7rs)/2, we search 
for the two of the range parameters, d^ N ^ N and c^-s^s, so as to reproduce 
the complex value of / = KN scattering length which was obtained by 



Martin's analysis; 0^(7=0) = —1-70 + z0.68 fm [26]. The range parameters 
with the calculated KN scattering amplitude are listed in Tabled] for f n = 90 
MeV case. For all combinations of kinematics and meson-baryon potential, 
we can find the range parameters which reproduce the Martin's value well. 
Also in cases of f n = 100, 110 and 120 MeV, we find such range parameters 
as given in Table IC.6I 

Using the meson-baryon potentials with these range parameters, we cal- 
culate the scattering amplitude in the 1 = channel. Figs. [3]l5]show the KN 
and 7rS scattering amplitudes for f n = 90 MeV. Fig. [3] shows an instructive 
result where the KSW potential is used in the NR kinematics. As seen in 
the figure, both the scattering amplitudes behave very singularly near the 
7rS threshold. We consider that such a singular behavior is caused by the 
mismatch of the kinematics and the potential. Since the flux factor in the 
KSW potential is relativistic, it should be used under relativistic kinematics, 
(at least, semi-relativistic one.) Indeed, such a singularity disappears in the 
cases (i)-(iii) where the potential is used with appropriate kinematics. (See 
Figs. H and [5]) 

As seen in Figs. H] and [5] the NR and SR give quantitatively the same 
KN scattering amplitude, though the 7rS scattering amplitude is different 
between two kinematics. Let us see the 1 = scattering amplitudes in more 
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Table 1: Range parameters for 7 = KN-ttT, system. f„ = 90 MeV case. d^ N g N 
{dns, tte) means the range parameter of Gaussian potential in KN-KN (7rE-7rS) channel. 
a KN (i=o) ^ s th° 7 = K N scattering length calculated with given range parameters. All 
quantities arc in unit of fm. 



Case 




(i) 


(ii) 


(iii) 


Kinematics 




Non-rela. 




Semi-rela. 


Potential 


KSW 


NRvl 


NRv2 


KSW 


d>KN, KN 


0.593 


0.576 


0.574 


0.487 




0.541 


0.725 


0.751 


0.457 


Re clk N ( 7=0 ) 


-1.701 


-1.700 


-1.700 


-1.703 


Im a^N (/=o) 


0.679 


0.677 


0.687 


0.677 
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Figure 3: 7 = scattering amplitudes calculated with the original KSW potential and NR 
kinematics. Left (right) panel shows KN (7rS) scattering amplitude. The real (imaginary) 
part of scattering amplitude is drawn with a black-solid (red-dashed) line. Vertical dashed 
line means the KN threshold. / ff = 90 MeV case. 
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Figure 4: I = scattering amplitudes for NRvl and NRv2 cases. The scattering ampli- 
tudes of NRvl (NRv2) are shown with thin (bold) line. Left (right) panel shows KN (7rE) 
scattering amplitude. /„. — 90 MeV case. 




| | ~ i l | 44 , 

1350 1400 1450 1350 1400 1450 

E [MeV] E [MeV] 



Figure 5: I = scattering amplitudes for SR case. Left (right) panel shows KN (7rE) 
scattering amplitude, — 90 MeV case. 
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Figure 6: 1 = scattering amplitudes for NRv2 case, calculated with various / w val- 
ues. Blue, green, orange and yellow lines correspond to f„—90, 100, 110 and 120 MeV, 
respectively. Left (right) panel shows KN (nT,) scattering amplitude. 



Table 2: Resonant energies estimated from the scattering amplitudes and the 7rS phase 
shift. The latter is estimated by 8^{E) = § and S'^(E) > 0. f„ = 90 MeV case. The 
energies are given in unit of MeV. 



Case 




(i) 


(ii) 


(iii) 


Kinematics 




Non-rela. 




Semi-rela. 


Potential 


KSW 


NRvl 


NRv2 


KSW 


E at Re f^ N = 


1402.5 


1408.4 


1411.1 


1417.0 


E at Re / wE = 


1354.6 


1388.8 


1393.0 


1422.0 


E at = ^ 


1354.6 


1388.8 


1393.0 


1422.0 



detail. In two NR cases (i) and (ii), the scattering amplitudes have the same 
behavior essentially. In the SR case (iii), the real part of the amplitude 
is smaller in the region far below the KN threshold, compared to the NR 
cases. We investigate the f n dependence of the scattering amplitude. Figs. 
|6] and [7] are the scattering amplitudes calculated with the NR and SR cases, 
respectively, when f n is varied from 90 to 120 MeV. In the NR case the 
scattering amplitude is found to depend strongly on the /„. value, especially 
near the 7rE threshold. It tends to be less attractive with larger value of f n 
used. On the other hand, in the SR case, it doesn't depend on the f w value 
so much. 

In all cases, the resonant structure is found in the KN and 7rE scattering 
amplitudes below the KN threshold. The resonant energies estimated from 
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Figure 7: I — scattering amplitudes for SR case, calculated with various f w values. Left 
(right) panel shows KN (ttT,) scattering amplitude. 



the scattering amplitude are summarized in Table HJ where they are consid- 
ered as the energy at which the real part of scattering amplitude becomes 
zero. The resonant positions on the KN scattering amplitude are around 
1410 MeV for NR case and near 1420 MeV for SR case. They are rather 
close to each other, (sightly lower in the NR case compared to the SR case.) 
On the other hand, on the ttT, scattering amplitude the resonant positions 
are rather different between two kinematics: about 1390 MeV in NR whereas 
about 1420 MeV in SR, for f w = 90 MeV. When the f w value is modified up 
to 120 MeV, the resonant position on the 7rE amplitude becomes higher to 
1417 MeV in the NR. Compared to the NR case, both the resonant positions 
on KN and ttT, amplitudes are stable for /„. value in the SR. (See Fig. |H]) 

3.2. Property of the 1 = KN-ttY, resonant state 

Using the meson-baryon potential determined in the previous section, we 
investigate the resonance in the 1 = KN-ttY, system. In practice, we search 
poles on the complex energy plane with the usual complex scaling method 
as explained in the section 12.21 

In all NR and SR cases, one pole is clearly foundjfl We denote the complex 
energy of a resonant pole as (E R , —T/2). The found poles for the case f n = 90 
MeV are shown in Table [3j The real part of energy E R is well determined to 



2 Preceding studies based on the chiral SU(3) theory reported that there exist two poles 
in the 1 = channel 17, 18| We will discuss a signature of another pole later. 
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Figure 8: f„ dependence of the resonance position appeared in scattering amplitudes. 
(Left) obtained from the KN scattering amplitude. (Right) obtained from the 7rE one. 
A symbol of black (blue) circle means NRvl (NRv2) and a symbol of red diamond means 
SR. 



Table 3: Pole position of the I — KN-irT, system and the meson-baryon distance in the 
pole state. {Er,—T/2) indicates the complex energy of the resonant pole. B^ N means 
the binding energy that is the Er measured from the KN threshold. The unit of these 
energies is MeV. (Vt^)kN' (v^ttE and (y/r*) kn+ttT, indicate the meson-baryon mean 
distance of KN, 7rS and total components, respectively. These values are given in unit of 
fm. = 90 MeV case. 



Case 


(i) 


(ii) 


(iii) 


Kinematics 




Non-rela. 


Semi-rela. 


Potential 


NRvl 


NRv2 


KSW 


Er 


1419.8 


1419.9 


1419.0 




15.2 


15.1 


16.0 


r/2 


26.0 


23.1 


14.4 




1.20- 0.35z 


1.25 - 0.35i 


1.21 - 0.49i 


vV>. s 


0.33 + 0.24i 


0.34 + 0.23i 


0.12 - 0.09z 


V / ( 7 ' 2 ) KN+tvT, 


1.20 - 0.28i 


1.25 - 0.29i 


1.21 - 0.49i 
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be about 1420 MeV, independently of kinematics and potential types. The 
imaginary part of energy T/2 depends on the kinematics; T/2 ~ 14 MeV in 
the SR, whereas T/2 ~ 25 MeV in the NR. Fig. [9] shows the pole positions 
of all cases when the /„- value is varied from 90 MeV to 120 MeV. In the NR 
cases, En is stable for /„- and is 1417-1420 MeV, but T/2 is ranging from 14 
MeV to 26 MeV. In the SR case the poles are found to distribute in compact 
region. The pole position (Er, —T/2) is determined with small deviation; 
(1419.5 ±1, —13 ±2) MeV. Thus, in the NR and SR kinematics the real part 
of energy of the pole is almost identical, though the imaginary one is different 
between two kinematics. The imaginary part of energy indicates a half width 
decaying to 7rS. In our study we have a constraint for KN channel but no 
constraint for ttT, channel. We consider that the difference of T/2 between 
two kinematics is due to the lack of constraint condition for 7rS channel. 

We have investigated another pole which is expected to exist, because 
many studies of the 1 = KN-nH system reported that this system has 
a double pole structure when an energy-dependent chiral SU(3) potential 



is used as we use |17J, |18j; The higher pole state is slightly below the KN 
threshold and with small width, while the lower one is far below the KN 
threshold and with large width, (for instance, the former is around (1432, 



-17) MeV and the latter is around (1400, -80) MeV |18|.) In our study, 
certainly we found self-consistent solutions which seem to indicate a lower 
pole of the double pole; the complex energy is (~ 1360, —40 ~ —90) MeV 
in case of NRv2. However, the separation of the pole from continuum states 
indicated by the 2#-line on the complex-energy plane seems insufficient and 
the pole position, in particular its imaginary energy, rather depends on the 
scaling angle 9. We consider that this is due to the limitation of the numerical 
accuracy. When poles have large decay width compared to excitation energy, 
it is known empirically that such poles are difficult to be found by the CSM 
with Gaussian base, since the spatial oscillation of the complex-scaled wave 
function is not well described with them. Therefore, we cannot conclude yet 
that the double pole structure is confirmed in the present our analysis. We 
need more investigation of the lower pole. 

We are interested in the inertial structure of A(1405). It has been in- 



vestigated theoretically in various ways [29j, [30|, |31|]. We consider the wave 
function of the pole state obtained by our ccCSM calculation. Fig. [10] shows 
complex-scaled wave functions of each component when the scaling angle 
9 = 30°. Here, the wave functions are multiplied by an appropriate phase 
factor so that the KN wave function becomes real at r = 0. Without this 
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Figure 9: I = pole in the complex-energy plane, calculated with NRvl, NRv2, and SR 
for various f n values. The number shown in the panel is the / w value which increases 
continuously from 90 to 120 along the connected dashed line in each case. 



phase factor, the complex-scaled wave function of KN-ttT, is normalized as 
Eq. (112p . Both of KN and 7rS wave functions are confirmed to be well 
localized. It is noted that localization of the 7rS component is due to the 
complex scaling, in spite that the state is above 7rS threshold. 

The mean distance between meson and baryon in the resonant state is 
calculated as 



MB 



) = (®e\r 2 MBt0 \$o), 



(34) 



where f MB = r Mc son ~ ^Baryon and $6/ means the complex-scaled wave func- 
tion of the resonant pole. It should be noted that the matrix elements of 
resonant state are obtained independently of 6 because the properties of 
the resonant wave functions are uniquely determined in the complex scaling 



method 32]. The expectation value calculated with Eq. (|34|) is inevitably a 
complex number because the resonant state is treated as a Gamow state in 
our framework. Therefore, the root-mean square distance, ( f MB ) , is also 
a complex number. Certainly its physical meaning is still unclear, but we 
show this quantity as a reference for the size. We believe that it is useful for 
us to get a feeling of the size of the system. Indeed, since the imaginary part 
of the obtained complex mean distance is smaller than its real part as will 
be shown later, we consider that the real part can be regarded as a mean 



distance with a physical meaning [32 
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Figure 10: Complex-scaled wave function of each component in the 1 = pole state at 
= 30°. We choose a phase such that the KN wave function becomes real at r — 0. Left: 
NRv2, Right: SR. /„ = 90 MeV case. 



If such an interpretation for the complex- valued distance calculated within 
the ccCSM is accepted, the mean distance between meson and baryon is 
about 1.2 fm in both cases of NR and SR kinematics with /„. = 90 MeV, as 
shown in Table EJ Those for other f„ cases are given in Table IC.81 In the SR 
case the mean distance has small f n dependence, remaining about 1.2 fm. In 
the NR case it depends on the /„. value similarly to the pole position, but it 
increases slightly to be 1.4 fm when f n is modified to 120 MeV. Thus, it is 
found that the mean distance is almost the same in both two kinematics. As 
for the imaginary part of the complex-valued distance, it is certainly small 
value of 0.5 fm at most. 

Compared to results of other studies, the size obtained by our calculation 
seems rather small, even if the modulus of a/ ( f\ 1B ) is regarded as a mean 
distance between meson and baryon in the 1 = system. For example, ac- 
cording to a study using a chiral SU(3)-based potential, it is 1.9 fm with 12 
MeV binding energy of KN (M = 1423 MeV) Q. We think that the differ- 
ence of size between two calculations is mainly caused by different definition 
of the resonant state. In the previous study, it is treated as a KN bound 
state, as a result of elimination of 7rE channel and perturbative treatment of 
the imaginary part of potential. On the other hand, the resonant state is a 
Gamow state in the current study since the complex scaling method imposes 
the outgoing boundary condition on a solution implicitly. 
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3.3. 1 = 1 KN-nE-nA system 

We make the same investigation on the isospin 1=1 sector which has 
three channels of KN, 7rS and 7rA. In our model of the Gaussian form 
potential, there are six range parameters in this sector. However, two of 
them, g^a^e and g^a^A) give no contribution to the result, since the potential 
strength of these channels are forced to be zero due to the SU(3) algebra that 
appears in the lowest-order potential of chiral SU(3) dynamics. (See C^ I=1 > 
in Eq. <Q) Since we assume d^N^T. = (d^ N j^ N + c? 7r2j7r2 )/2 similarly to the 
1 = case, three range parameters, d^ N j^ N , d^^s and d^N^A, are unknown 
parameters to be determined. 

Similarly to the 1 = case, we constrain the range parameters in the po- 
tential by the Martin's value of the 1 = 1 KN scattering length; cikn(i=i) = 



0.37 -MO. 60 fm 26]. However, the three unknown parameters can't be deter- 
mined by only the complex value of (1kn(i=i)- Here, we reduce the number 
of unknown parameters by referring the following two facts: 1. In studies 
with chiral unitary model, isospin symmetric subtraction constants have been 
often assumed and succeeded to reproduce various physical quantities [33 . 
2. In a separable potential used in Faddeev-AGS calculation of KNN-irYN 
j5[, the cut off parameter for the KN channel is not so different between 
1 = and 1=1 sector. Based on these facts, we examine three conditions 
as follows: 

- Cond. (a) d^N^RN is fixed to that of the 1 = case. 

^ttE.ttE and d^N^A are searched to reproduce the complex 
value of o>kn(i=i)- 

- Cond. (b) d^N,KN an d g^e^e are fixed to those of the 1 = case. 

dRN,wA is searched to reproduce the real part of clkn{i=x)- 

- Cond. (c) Similar to the condition (b), but d^N^A is searched to 

reproduce the imaginary part of a^ N ^ I=1 y 
We describe mainly the result with f„ = 90 MeV. In the NR case, we 
can find a set of range parameters which satisfy each condition (a)-(c), as 
shown in Table HI However, it is found that the 1 = 1 scattering amplitude 
calculated with the condition (a) has a resonant structure slightly below 7rS 
threshold as shown in the left column of Fig. [Til although no resonance 
in / = 1 sector has been confirmed theoretically and experimentally. In the 



3 In this case we have found a pole on the complex-energy plane by the ccCSM at 
(E R ,-T/2) = (1326.4,-1.3) MeV. Certainly, this pole exists only by 5 MeV below ttE 
threshold (=1331 MeV). 
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Table 4: Range parameters for the 1 = 1 KN-tt'S-ttA system with the non-relativistic 
kinematics. f„ = 90 MeV case. The same table as that in Table [TJ "Condition" is 
explained in the text. 



Case 




(i) 






(ii) 




Kinematics 
Potential 




Non-rela. 
NRvl 




Non-rela. 
NRv2 




Condition 


(a) 


(b) 


(c) 


(a) 


(b) 


(c) 


URN, KN 
dRN, ttA 


0.576 
0.212 
0.355 


0.576 
0.725 
0.502 


0.576 
0.725 
1.038 


0.574 
0.221 
0.391 


0.574 
0.751 
0.558 


0.574 
0.751 
1.138 


Re (Xk N ( /=1 ) 
Im Or-jv(j=i) 


0.372 
0.601 


0.371 
1.444 


0.650 
0.600 


0.371 
0.601 


0.370 
1.358 


0.621 
0.600 



conditions (b) and (c), such a resonant structure does not appear in all the 
scattering amplitudes. (See the right column of Fig. [TT]) 

In the SR set of range parameters to satisfy the conditions is found 

only in the case (a), but is not found in the cases (b) and (c). (See the left 
three columns in Table |5]) Scattering amplitudes calculated with the range 
parameters of the condition (a) is shown in the left column of Fig. [12] As 
seen in this figure, the 7rS scattering amplitude indicates repulsive nature of 
7rS-7rS channel, in spite that the direct 7rS potential is originally attractive. 
This is considered to be a consequence of the coupled- channel effect. 

We investigate other f n values such as 100, 110 and 120 MeV. Range 
parameters and KN scattering length for these f n 's are listed in Table 1(191 
It is confirmed that the scattering amplitudes for these f^s are essentially 
the same as those for f n = 90 MeV case above mentioned. 

We consider the cases where the conditions are slightly relaxed. First, 
we vary the value of the range parameter d^N,RN slightly in the condition 
(a), since it may be too strict constraint that d^ N j^ N is fixed to that of 
the 1 = sector. With 10% modification of the dRN,RN: essential property 
of the scattering amplitude is found to be unchanged; In the SR the 7rS 
scattering amplitude still indicates the repulsive nature, and in the NR a 
resonant structure is kept to appear around the 7rS threshold. Next, we 
consider that only the imaginary part of (3^(7=1) is reproduced in 

the condition (a) with the SR case, giving up the reproduction of its real 
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Figure 11: 1 = 1 scattering amplitude calculated with NRv2. Left (Right) panels are 
calculated with the condition (a) (condition (c)). Top: KN — > KN, Middle: 7rE — > 7rE, 
Bottom: 7rA -> ttA. = 90 MeV case. 
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Figure 12: 1=1 scattering amplitude calculated with SR. Left (Right) panels are calcu- 
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Bottom: 7rA ttA. f n = 90 MeV case. 
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Table 5: Range parameters for the 1=1 KN-tt'E-ttA system with the semi-relativistic 
kinematics. /„■ = 90 MeV case. Columns "(a)-Rl" to "(a)-R3" show the results obtained 
with the condition (a) relaxed. Details are explained in the text. 



Case 




(iii) 




(iii) 




Kinematics 
Potential 




Semi-rela. 
KSW 




Semi-rela. 
KSW 




Condition 


(a) 


(b) (c) 


(a)-Rl 


(a)-R2 


(a)-R3 


d>KN, KN 

dRN, tv a 


0.487 
0.313 
0.422 


0.487 0.487 
0.457 0.457 
0.187 0.100 


0.487 
0.960 
1.550 


0.487 
1.210 
1.000 


0.487 
1.340 
0.930 


Re a KN (i=i) 
Im axjv(/=i) 


0.374 
0.602 


0.072 -0.145 
0.268 0.404 


0.932 
0.591 


0.871 
0.600 


0.869 
0.600 



part. We can obtain some sets of range parameters satisfying such a relaxed 
condition (a) as listed on the right three columns in Table |5j When these 
range parameters are used, the 7rS scattering amplitude is found to indicate 
attractive nature as seen in the right column of Fig. [T2] Compared with the 
NR result shown in the right column of Fig. [TT] (NRv2-c), it is noticed that 
the SR and NR give quite a similar result for all channels of KN, 7rE and 
7rA, in such that Im n (i=i) is fixed to be the experimental value. 



4. Summary and future plan 

We have studied a KN-ttY system with a coupled-channel complex scal- 
ing method (ccCSM) [23J using a chiral SU(3) potential. In our study, 
scattering states as well as resonant states are investigated within a sin- 
gle framework of ccCSM. Resonant poles are obtained by diagonalizing a 
complex-scaled Hamiltonian with Gaussian base, similarly to bound states 
calculation. Scatte ring problem is solved with an advanced use of ccCSM, 



ring 

"CS-WF" method 25j. In the CS-WF, due to Cauchy's theorem scattering 
amplitudes are calculated with a complex-scaled wave function which is also 
described with Gaussian base. Thus, both of resonant and scattering prob- 
lems can be solved with Gaussian base and therefore they can be treated with 
small modification of bound-state calculation. This is the most advantageous 
point of ccCSM. 
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Based on Ref. 14J where a meson-baryon potential is derived from a 



chiral SU(3) theory, we have constructed a meson-baryon potential (KSW 
potential) which is a local potential with Gaussian form in r-space. Since in 
the present study it is necessary to deal with pion whose mass is very light, 
we have examined semi-relativistic kinematics (SR) as well as non-relativistic 
one (NR). In the NR case, non-relativistically approximated versions of the 
KSW potential are used. 

By using the CS-WF method, range parameters of our Gaussian-form 
potential are determined for both the NR and SR kinematics, so as to repro- 
duce the value of KN scattering length obtained by Martin's analysis [26] . 
The scattering amplitudes are investigated with the CS-WF method using 
the determined potentials. 

In the 1 = sector, it is found that the KN scattering amplitude near the 
KN threshold is well constrained by the KN scattering length, since its f n 
dependence is small and both kinematics give similar results. However, far 
below the KN threshold the scattering amplitudes are strongly dependent 
on the fx value in NR and the results of NR and SR are rather different from 
each other. We consider that further data far below KN threshold, such as 
7rS scattering length, are necessary to reduce such an uncertainty in the deep 



KN bound region, as pointed out in Ref. |34j. A resonant structure is seen 
below the KN threshold in the 1 = KN and 7rS scattering amplitudes. 
For the modification of the f n value, it is confirmed that in the NR case the 
resonant position of 7rS scattering amplitude is rather varied, whereas in the 
SR case it is quite stable around 1420 MeV in both of KN and 7rE scattering 
amplitudes. 

Similarly, the 1 = 1 sector has also been investigated. When the potential 
of our model is constrained by the complex value of I = 1 KN scattering 
length, it is found that the scattering amplitude has a resonant structure 
slightly below the 7rE threshold in the NR case, and that 7rS scattering 
amplitude shows repulsive nature in the SR case. However, if the constraint 
condition for the potential is relaxed so that only the imaginary part of the 
KN scattering length is reproduced, such a resonant structure disappears in 
the NR case and the 7rS scattering amplitude becomes attractive in the SR 
case. 

Property of the resonant pole in the 1 = sector corresponding to the 
A(1405) has been studied with the usual ccCSM. For f n = 90 ~ 120 MeV, 
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resonance pole is found around 



(E R , -r/2) = { ( 




on the complex energy plane. In the NR case, imaginary part of the pole 
energy is rather dependent on the f n value, while in the SR case, both real 
and imaginary parts of the pole energy are quite stable for the value. 
As well as the pole energy, it is found that the f n dependence of various 
quantities is small in the SR calculation rather than in the NR one. We have 
estimated the "size" of the pole state by evaluating a root-mean distance with 
a bi-orthogonal set of complex-scaled wave function, \/ ( f 2 MB ) . This quantity 
does not mean the physical meson-baryon distance since it is natural to be 
a complex value, but is expected to give us a guide of mean distance. The 
calculated meson-baryon "distance" s are similar value for the NR and SR 
cases; 1.3 — «0.3 fm for the NR case and 1.2 — z0.5 fm for the SR case. 

As a result of the current study with the NR and SR kinematics, it is found 
that both the kinematics give essentially the same result on KN quantities 
such as the KN scattering amplitude, the real energy of / = pole state 
and its size, when we constrain the model by the KN scattering length that 
is the quantity at the threshold. However, difference between them becomes 
large as for the quantities far below the KN threshold and those related to 
7rS, where the relativistic effect becomes important. 

The pole above discussed is considered to be the higher pole of the double 
pole of A (1405). As mentioned in the section 13.21 we found a signature 
of the lower pole around (~ 1360, —40 ~ —90) MeV in a non-relativistic 
case. However, we can't conclude that this is the lower pole, because its 9 
trajectory is somehow unstable. We think that this is due to limitation of 
numerical accuracy of the CSM with Gaussian base for the case of large decay 
width. The poles of broad resonances can be investigated by applying an 
analytic continuation in the coupling constant to the complex scaling method 
(ACCC+CSM) [35]. It is one of our future plans to carry out ACCC+CSM 
and clarify whether our potential leads to the double-pole structure or not. 

Thus, we have a KN-nY potential for both isospin channels, which is 
based on a chiral SU(3) theory and is a local Gaussian form in r-space. In our 
future plan, we will investigate the three-body system of K~pp (KNN-7iYN 
system with J 71 " = 0~, T = 1/2) which is the most essential kaonic nuclei. 
Since the ccCSM can adequately deal with resonant states of a multi-channel 
system in principle and the CSM is known to be effective for the nuclear 
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many-body study 24| , we expect that the ccCSM will give a pole position of 
the K~pp accurately and reveal its structure. We have interested in the role 
of tiYN three-body dynamics, because its implicit/explicit treatment may 



cause a large difference in the binding of K~pp as pointed out in Ref. (36 
It is expected that the contribution of nYN three-body dynamics will be 
investigated with the ccCSM. 

We think that it is worthwhile to use the updated value of the KN scat- 
tering length, instead of the Martin's value obtained from old data, in our 
analysis. The SIDDHARTA group reported the K~p scattering length di- 
rectly from the precise measurement of the shift and width of the Is atomic 



level energy of kaonic hydrogen atom [20] . Using this data, the K n scatter- 



ing length is estimated with the coupled-channel chiral dynamics 22] . These 



values of K~p and K~n scattering lengths are available in our calculation. 
Furthermore, the SIDDHARTA group is planning to perform experiments 
on kaonic deuterium in the SIDDHARTA-2 experiment. These forthcoming 
experiment will provide us the more precise values of KN scattering length 
for both isospin channels 37 . 

Furthermore, the ccCSM approach can be applied to other hadronic sys- 
tems. For instance, it seems interesting to investigate with this method a 
few-body system involving D meson in the charm sector, which is an analo- 
gous system with K meson in the strangeness sector [38| . 
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Appendix A. Scattering amplitude for a multi-channel case 

We consider a Schrodinger equation for a multi-channel system; H\ty) = 
E\^f). Hamiltonian H and total wave function are given as 

H = $> c °|c>(c| + £ Vcc>\c)(c'\ } |*) = l*c>|c>, (A.l) 

c c,c' c 

where H® is the kinetic energy operator p 2 c /2[i c for the channel c and V cc r is 
a potential between channels c and c'. |\I/ C ) is a wave function of a channel 
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c. Projecting the Schrodinger equation onto a channel c, we get such an 
equation for the c-channel wave function as 

H°\V C ) +Y,V CC <\V C <) =E\* C ). (A.2) 

c,c' 

When the incident channel is cq and the incoming wave function is given as 
|0c o ,fc co ); the above equation can be modified formally as 

l^ co) > = \<Pc lhkco )S cco + E _^ + ie J2 V ^ ) ^ (A.3) 

c d 

taking into account the outgoing boundary condition appropriately. Here, we 
write the incoming channel Cq on the wave function of each channel explicitly. 
The above equation is expressed in r-space as 

^ CO) (^)=0 Co , fcco (r Co )5 CCo 

+ E [ d < G e (r e y c ;£0 (r' c \V cc ,\¥^), (A.4) 

where G c (r c , r' c ; E) is a Green function in the channel c. By performing 
complex integral as shown in many textbooks, it becomes 

1 1 Oil p ik c\ r c- r 'c\ 



c 



1 9 II P lfe cr c 

_/ ' (\r c \ » \r>\), (A.6) 



An h 2 r c 

where k c = k c r c /r c and r c = |r c |. Then, the channel-c wave function becomes 
* ( c°\r c ) = MrJ5 CC0 + ^ x (-±-^j ^^(^JKc'I^}, (A.7) 

using the fact that the function Ci fc c (r c ) is e lkc ' Vc . Substituting the channel-c 
wave function in Eq. I 1A.1I) with the above expression, the total wave function 
is given as 

l^> = ^c (^ co )|co> + ^ ^|c> ■ (-^) ^ 5Z(^, fc J^I^ 0) >- (A.8) 

c c ^ ' d 
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Therefore, the scattering amplitude between the initial channel Cq and the 
final channel c is 

f CCo (k c ,k Co ) = (A.9) 

c' 

The wave functions <p c ,k c (, r ) = e lkc ' r and ^J°\r) are expanded on partial 
waves / as 

Clfcc (r) = ^J2 ilJj ^ Y i*rn(^)Y lm (n r ), (A.10) 

K c r 

*0*)(r) = 47r^i J ^^y4(Q feco )Y im (Q I .) (A.ll) 
and the scattering amplitude is expanded on the orbital angular momentum 



as 



f CC0 {k c ,k CQ ) = ^2(2l + l)P l (cose ikc>kco) )f ljCCO (k c ,k co ). (A.12) 
i 

Thus, expanding Eq. (1A.9|) for the partial waves using Eqs. flA.10D - fjA.12p . 
the scattering amplitude for the partial wave I is given as 

f, cco (k c ,k co ) = --^^^OiMlV^^ir)), (A.13) 



»c»oo c , 



if the potential V^ c / is a central potential. 

Appendix B. Matrix element of the kinetic term in the semi- 
relativistic case 

In the semi-relativistic kinematics, the kinetic energy and mass terms in 
the Hamiltonian are of the form of \J m? + p 2 + \J M 2 + p 2 as shown in Eq. 
(J2J). In this article, a wave function is expanded in terms of partial waves; 
^( r ) = J2im r_1 V'«( r )^m(^)- Furthermore its radial part is expanded with a 
Gaussian base, as Eq. (IH)l and Eq.( |24j) explained in non-relativistic case; 

Mr)/r = X) c 3 G 5( r )/ r > G K r ) = ^(&i)^ +1 exp[-r 2 /26f], (B.l) 
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where Ni(bj) means a normalization factor. We need to calculate the matrix 
element (r _1 Cr*Yz m | \Jm 2 + p 2 \r~ 1 G l jYi/ m /). This matrix element is calculated 



as follows: 



l -G\Y lm Wmi + pi \^G l -Y Vm ,) 
drdr'dqdq' l^G\Y lm \r) (r\q) (q\V™ 2 + p*\q') (q'\r') (r'^-G*) 

POO 

/ dqq 2 \/ m 2 + q 2 
Jo 



Su'5. 



ll'Umm' 



x / drG^jtiqr) x / dr'G)(r')ji{qr'), 



(B.2) 



where (r\q) = e tqr and its expansion (Eq. (lA.lOj) ) is used. In case of s-wave 
(/ = 0) which we are considering in this article, the last integration for r and 
r' variables can be performed analytically. Finally, the above matrix element 
for a complex-scaled p is expressed as 



(^YooWm 2 + (pe^) 2 | ±G°Y 00 ) 

-^(&A) 3/2 / dqqWm 2 + q 2 e- 2it> exp 
V 71 " Jo 



(B.3) 



The integration for the variable q is carried out numerically. The matrix 
element of Eq. (1B.3|) is used when we diagonalize the complex-scaled Hamil- 
tonian to find resonant poles as explained in section 12.21 and we calculate 
a complex-scaled wave function as shown in Eq. ( 1251) to obtain scattering 
amplitudes. 

Note that in a semi-relativistic case any wave number k is calculated 
from \j m 2 + hk 2 + a/ M 2 + hk 2 = E and that any reduced mass /i in a 
non-relativistic case is replaced with corresponding reduced energy u = 
E M E B /(E M + E B ) where E M and E B mean meson and baryon energies, 

respectively; Em = V m 2 + Kk 2 and E B = \/M 2 + hk 2 . 



Appendix C. Detailed results of f n = 100, 110, 120 MeV cases 

We show the results for the cases of /„- = 100, 110 and 120 MeV. Tables 
IC.6t IC.7I and IC.8I correspond to Tables [fl |2] and |3] which are for the case of 
/„- = 90 MeV, respectively. Table IC.9I shows the result corresponding to the 
case (ii) in Table H] and the cases (iii)-(a) to (c) in Table [5j 
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Table C.6: Range parameters for / = KN-ttT, system. All quantities are in unit of fm. 



f v = 100, HQ and 120 MeV cases. Corresponding to Table [J 



Case 


(i) 


(ii) 


(iii) 






Kinematics 




Non-rela. 


Semi-rela. 






Potential 


NRvl 


NRv2 


KSW 


/* = 


100 














URN, KN 


0.503 


0.501 


0.421 








0.665 


0.695 


0.395 






Re a_R-iv(i=o) 


-1.702 


-1.702 


-1.703 






Im a KN (i=o) 


0.680 


0.683 


0.673 


f* = 


110 














d<KN, KN 


0.440 


0.438 


0.369 






«7r£,7r£ 


0.605 


0.636 


0.348 






Re a^7v(/=o) 


-1.701 


-1.700 


-1.696 






Im a^ N ( /=0) 


0.681 


0.681 


0.681 


u = 


120 














">KN, KN 


0.386 


0.384 


0.327 






«7rE,7rS 


0.549 


0.581 


0.310 






Re a^7v(/=o) 


-1.699 


-1.700 


-1.690 






Im clr n ( /=0 ) 


0.674 


0.669 


0.669 
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Case 




(i) 




(ii) 




(iii) 






Kinematics 






Non-rela. 




Semi-rela. 






Potential 




ATT~) 1 

JNKvl 




ATT) O 




KSW 


u = 


100 


Er 

Brn 
1 /A 




1417.3 
17.7 
23.1 




1418.0 
17.0 
19.8 




1419.6 

_L<J.4: 

13.2 








l.ol 


— U.zyz 


l.OO 


— U.oOt 


1.21 


- 0.49i 


f* = 


110 


Er 

Ern 
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16.6 
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Table C.9: Range parameters for / = 1 KN-tt'S-ttA system with non-relativistic kinematics 
(left three columns) and semi-relativistic kinematics (right three columns). / w = 100-120 
MeV case. Corresponding to Tables [4] and [5] 



Case 




(ii) 






/ • • • \ 

(m) 




Kinematics 




Non-rela. 




Semi-rela. 




Potential 






NRv2 




KSW 




Condition 


(a) 


(b) 


(c) 


(a) 


(b) 


(c) 


fn = 100 
















URN, KN 




0.501 
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0.421 


0.421 


0.421 
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0.695 


0.277 


0.395 


0.395 
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0.348 


URN, ttA 




0.221 
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0.372 
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Im a,R N (/= 


=D 
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1.504 
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0.600 


0.154 


0.337 


A = 120 




















0.384 


0.384 


0.384 


0.327 


0.327 


0.327 






0.139 


0.581 


0.581 


0.223 


0.310 


0.310 


«/?iV, TTA 




0.173 


0.240 


0.324 


0.236 


0.140 


0.105 


Re a^^y (/= 


=i) 


0.382 


0.369 


0.689 


0.379 


0.123 


-0.097 


Im cikn (i= 




0.608 


1.539 


0.601 


0.600 


0.141 


0.360 
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